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Postprocessing can speed up general quantum search algorithms 
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Department of Physics, IIT Bombay, Mumbai-400076, Indi£| 

A general quantum search algorithm aims to evolve a quantum system from a known source state 
|s) to an unknown target state \t). It uses a diffusion operator Da having source state as one of its 
eigenstates and It, where denotes the selective phase inversion of \ip) state. It evolves |s) to a 
particular state \w), call it w-state, in 0{B/a) time steps where a is |(t|s)| and B is a characteristic 
of the diffusion operator. Measuring the w-state gives the target state with the success probability 
of 0{ljB^) and O(B^) applications of the algorithm can boost it from 0(1/B^) to 0(1), making 
the total time complexity 0{B^/a). In the special case of Grover’s algorithm, Dg is 7s and B is very 
close to 1. A more efficient way to boost the success probability is quantum amplitude amplification 
provided we can efficiently implement Such an efficient implementation is not known so far. In 
this paper, we present an efficient algorithm to approximate selective phase inversions of the unknown 
eigenstates of an operator using phase estimation algorithm. This algorithm is used to efficiently 
approximate Iw which reduces the time complexity of general algorithm to 0(B/a). Though 0{B/a) 
algorithms are known to exist, our algorithm offers physical implementation advantages. 

PACS numbers: 03.67.Ac 


I. INTRODUCTION 


Grover’s algorithm or more generally quantum am¬ 
plitude amplification drives a quantum computer from 
a source state |s) to a target state \t) by using selec¬ 
tive phase inversion operators. Is and It, of these two 
states [in. The algorithm is an 0{l/a) times itera¬ 
tion of the operator A{s,t) = —Iglt on |s) to get |t) 
where a = |(t|s)|. For search problem, we choose |s) 
to be the uniform superposition of all N basis states to 
be searched i.e. |s) = \i)/^fN. In case of a unique 

solution, the target state |7) is a unique basis state and 
a = |(t|s)| = Thus Grover’s algorithm outputs 

a solution in just O^s/ll) time steps which is quadrati- 
cally faster than classical search algorithms taking 0{N) 
time steps. Grover’s algorithm is proved to be strictly 
optimal Q. 

A general framework of quantum search algorithms 
was presented in Q. The general quantum search al¬ 
gorithm replaces Is by a more general diffusion operator 
Ds with the only restriction of having |s) as one of its 
eigenstates. This restriction seems to be more or less 
justified as the diffusion operator should have some spe¬ 
cial connection with the source state. Let the normalized 
eigenspectrum of Ds be given by Ds\i) = with 

\t) as the eigenstates and {Of) as the corresponding 
eigenvalues (eigenphases). Since a global phase is irrel¬ 
evant in quantum dynamics, we choose 77^,1 s) = |s), i.e. 
Of^s = 0. The general search operator is 5 = Dsit and 
its iteration on |s) can be analyzed by finding its eigen¬ 
spectrum. Such an analysis was done in Q and we found 
that the performance of general algorithm depends upon 
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two quantities Ai and A 2 , where 

= ( 1 ) 

is the p*** moment of cot with respect to the distri¬ 
bution |(7|t)p over all 7 s. For a successful quantum 
search, we found Ai = 0{aB) « 0 to be an essential 
condition. Here B = \/l + K 2 is a characteristic of the 
diffusion operator Ds- An algorithm was also presented 
in Section IV.A of @ which uses an ancilla qubit to clev¬ 
erly control the applications of Ds and D\ to effectively 
make Ai zero. In this paper, we restrict ourselves to the 
case Ai = 0. In case, it is not so, we can always use 
just-mentioned algorithm to make it so. 

Unlike Grover’s algorithm, the general algorithm does 
not rotate the quantum state in a plane spanned by |s) 
and \t). Rather, it induces rotation in a different plane 
spanned by |s) and a particular state \w), which we refer 
as w-state here. The rotation angle is 2a/B and w-state 
has an overlap of 0(1/77) with the desired final state \t). 
Thus, we get w-state after one round of algorithm which 
takes 0{B/a) time steps and we perform 0{B^) rounds 
to get \t) with high probability. That makes the total 
time complexity to be 0{B/a) x 0{B'^) = 0{B^/a). 

We know that performing 0{B^) rounds of algorithm 
to boost the success probability from 0(1/77^) to 0(1) 
is a classical process. Quantum mechanically, it can be 
done more efficiently using quantum amplitude amplifi¬ 
cation (QAA). But to drive [w) towards \t), we need to 
implement I^j, the selective phase-inversion of w-state. 
So far, Iw is not known to be efficiently implementable. 
By definition |w) = 5^|s) for some q so 1^ = S'^IsS~'^. 
So we can implement 1^ if we can implement Ig ■ But the 
whole purpose of studying general quantum search algo¬ 
rithm was to discuss the cases when Is is not efficiently 
implementable. Rather, what is available to us is a gen¬ 
eral diffusion operator Dg of which Is is just a special 
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case. We point out that Is is not easily implementable 
in cases of physical interest il- 

Thus, we are in some sense forced to use the classical 
process of running 0{B‘^) rounds of search algorithm to 
get the target state. However, in this paper, we show 
that Iw can be efficiently implemented and 0{B/a) time 
steps are enough to get the target state. We use quan¬ 
tum fourier transform to implement I^- We point out 
that 0{B/a) algorithms are known to exist do but 
our alternative approach offers implementation advan¬ 
tages as in general, it uses significantly lesser number of 
controlled transformations which are physically harder 
to implement. In next section, we present the algorithm. 
The crux of our algorithm is the approximation of Iw 
using phase estimation algorithm. We present the algo¬ 
rithm for this approximation in Section III. We conclude 
in Section IV. 


II. ALGORITHM 

We briefly review the general quantum search algo¬ 
rithm 0. It iterates the operator Dglf on |s). Here Ds 
is as defined earlier and if is the selective phase rotation 
of target state by angle of (j). We choose (/> = tt so that 
if = It- We assume |s) to be a non-degenerate eigenstate 
of Ds for simplicity. The normalized eigenspectrum of Ds 
is given by Ds\€) = By convention, 9i=s = 0. Let 

other eigenvalues satisfy 

II > 6*111111 > 0, 6 i^G[- 7 r, 7 r] (2) 

The iteration of 5 = Dglt on |s) can be analysed by find¬ 
ing its eigenspectrum. The secular equation was found 
in @ to be 

^|(^|i)|2cot^^^ = 0. (3) 

i 

Any eigenvalue of S has to satisfy above equation. 

Since cot x varies monotonically with x except for the 
jump from — oo to oo when x crosses zero, there is a 
unique solution Afc between each pair of consecutive 0^’s. 
As = 0, there can be at most two solutions Afc in 
the interval [—0min, ^min]- Let these two solutions be A±. 
We have |A±| < Smin- The eigenstates |A±) correspond¬ 
ing to these two eigenvalues are the only relevant 
eigenstates for our algorithm. Precisely, if we assume 
|A±| <C 0inin, then the initial state |s) is almost com¬ 
pletely spanned by |A±). The eigenphases A± are given 
by 

A± = ±^(tan77)=^i ; cot277 = . (4) 

where 

H = vTT^, Ap = ^|(€|t)|2cot^’| . (5) 

l^s 


We consider the case when Ai = 0. Then eq. o indi¬ 
cates that 

TT /fy 

Ai = 0^r;=-, A± = ±- . (6) 

With 77 = 71/4 and (/ = tt, Eq. (23) and (24) of Q 
gives us the initial state |s) and the effect of iterating S 
on |s) in terms of two relevant eigenstates |A±). We have 

|s) = -7/^[e*^+/^|A+)-e*^-/2|A_)], (7) 

and 

5«|s) =-7/V2[e*«'^+|A+)-e"«'^-|A_)], (8) 

where q' = q + \ . 

For q = q^a ~ 7r/2|A±| = 7rH/4a, the state iS'^“|s) is 
almost the w-state \w) given by 

S‘>-\s) = \w) = l/V2i\X+) + \X.)). (9) 

As shown in we have Ktlw)! = 1/B. More explicitly, 

|t) = |H + |A^), (10) 

where |Aj_) denotes any state orthogonal to |A±). If 
we can efficiently implement Iw then 0{B) iterations 
of quantum amplitude amplification operator A{w, t) = 
Iwh on jw) will bring us close to the target state \t). 

We make a simple but useful observation that for our 
purpose Iw is completely equivalent to /a± , the selective 
phase inversion of the two-dimensional subspace spanned 
by |A±). For our purpose, the quantum state is always of 
the form A'^{w,t)\w) for an integer q. As the amplitude 
amplification operator A{w, t) rotates the state vector in 
the plane spanned by \w) and |t), the state A'^{w, t)\w) is 
always a linear combination of jic) and \t). Now, consider 
the two mutually orthogonal basis states |±) = |A+) ± 
|A_) of |A±) subspace. We note that |-|-) = \w) and 
Eq. (fTUl) implies that \t) also have zero component in |—) 
state. Thus A‘^{w,t)\w) has zero component in |—) state. 
For such a state, the operator which is a selective 
phase inversion of both states |-|-) and |—), is completely 
equivalent to Iw, the selective phase inversion of only 
|-|-) = |w) state. Hence /a± is as good as Iw for using 
quantum amplitude amplification to evolve Iw) to \t). 

The operator is implemented using the phase esti¬ 
mation algorithm and our partial knowledge |A±| < 0min 
of the corresponding eigenphases. This knowledge helps 
us to distinguish |A±) from other eigenstates |A_l) of S. 
Next section shows that 0{\nB /queries are suffi¬ 
cient to approximate I\^ with an error of o{l/B). Here 
o(x) <C a; is the small-o notation used in computer science 
literature. As we need 0{B) applications of Iw = Ix±, 
the total error in our algorithm is Bo{l/B) -C 1 as er¬ 
rors accumulate linearly in quantum mechanics. The 
total number of queries required to evolve |w) to \t) is 
0{B) X 0{\nB/e^in) = O(BlnH/0„i„). 
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The total time complexity of our algorithm is the sum 
of (i) the time complexity 0{B/a) of first evolution from 
|s) to Ire) and (ii) the time complexity 0{BhiB/Orai-a) of 
second evolution from |w) to \t). This is 


O 



B\nB 

^min 


( 11 ) 


For typical cases when a <C 0min/ln-B, above term is 
0{B/a). 

Our algorithm is not the first 0{B/a) algorithm for 
general quantum searching but it offers implementation 
advantages over the earlier algorithms. To understand it, 
we consider the controlled quantum search algorithm pre¬ 
sented in @ having a time complexity of 0{B/a). But 
that algorithm succeeds by using 0{B/a) controlled ap¬ 
plications of the operators Ds and It ■ Similarly, recently 
a different 0{B/a) algorithm was presented in [T^ but 
that also uses 0{B/a) controlled applications of Dg as 
well as 0{l/a) controlled applications of It- In com¬ 
parison to these two algorithm, the alternative presented 
here uses only 0(i? In controlled applications of 

S = Dglt- These controlled applications are required for 
the Phase Estimation algorithm as we illustrate in the 
next section. For the typical cases, lni?/0min <C 1/a 
and hence the present algorithm uses significantly lesser 
number of controlled operations compared to earlier al¬ 
gorithms. As controlled operations are harder to imple¬ 
ment physically, our algorithm offers implementation ad¬ 
vantages. 


III. APPROXIMATE SELECTIVE PHASE 
INVERSION 

A. Basic Idea 

We present an algorithm to approximately implement 
the operator We attach a workspace with the 

Hilbert space Bw to our main quantum system on which 
we wish to implement We denote the Hilbert 

space of main quantum system by Hm and refer to it as 
mainspace. We work in the joint space Hj — Hm ® Hw ■ 
Let the initial state of this joint space be 

\(,^) = \0W), ( 12 ) 

where \a) is any standard known state of the workspace. 
We expand the initial state of mainspace |^) in the eigen- 
basis of the operator S as 

10 = (13) 

± h 

where |A±) are those two eigenstates on which we wish to 
implement the selective phase inversion and | denote 

any of the remaining eigenstates orthogonal to |A±). The 

action of operator I\j_ on the mainspace is given by 

(/a± ® 1™) lOo-) = I0)k), (14) 


where lu, is the identity operator acting on the workspace 
and 

10) = h±\0 = X! (“f±) l-^±) + '^Ch\X±,h)- (15) 

it h 

We point out that the index h is not necessary in our 
analysis and we can represent all |Aj_^;i) by |A_l) for con¬ 
venience. 

To approximate the operator I\^ on the mainspace, 
we aim to obtain the state |^'). For this, we partition 
the workspace Hw into two mutually complementary sub¬ 
spaces: Z and which can be easily distinguished from 
each other. So we have Hw = Z\J Z^. Suppose we can 
design an operator C whose action on the joint space is 
given by 

C(|A±)|a)) = |A±)|Z) + |0(c)), 

C(|A^)|a)) = |A^)|Z^) + |0(c)). (16) 

Here \Z) and \Z^) denote a state completely within the 
subspaces Z and Z^ respectively. Also, e <C 1 and |0(e)) 
denotes any arbitrary state of small length 0(e). Con¬ 
sider the following operator 

n = c^ {imZ>iz)c (17) 

where Iz is the selective phase inversion of .Z-subspace 
acting on the workspace and Im is the identity operator 
acting on the mainspace. As .Z-subspace is easily distin¬ 
guishable from its complementary subspace Z^, imple¬ 
mentation of is straightforward. 

Using Eq. (ITBl) . we find that 

{lm®Iz)C{\\±)\a)) = -|A±)|Z) + |0(e)), 
{lm^Iz)C{\Xj_)\a)) = |Ai)|Z^) + |0(e)). (18) 

Also, left multiplication of Eq. (HU) by cl and the uni- 
tarity of C gives us 

Ct(|A±)|Z)) = |A±)|a) + |0(e)), 
Ct(|A^)|Z^)) = |A^)|a) + |0(e)) (19) 

Eqs. (fT71) . (fT51) and (ITOl) imply that 

7e(|A±)|a)) = -|A±)|a) + |0(e)) 

7e(|Ax)k)) = |Ax)k) + |0(e)) (20) 

Using the linearity of TZ and Eqs. (I13I15I) along with the 
above equation, it is easy to check that 

7^ie,^) = ir)k) + io(e)) (21) 

Thus we can approximate the operator and 

hence upto an error of 0{e) provided we can design 
an operator C satisfying Eq. (HH). The tolerable value 
of e depends upon the total number of approximate im¬ 
plementations of required by the algorithm. In our 
case, this number is 0{B) and 0{B) approximate imple¬ 
mentations will contribute an error of 0{Be) as errors 
add up linearly. Thus for e <C 1/H, our approximations 
will contribute a negligible error. 
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B. Phase Estimation Algorithm 


where the state \(j)\) can be written using Eq. (l26t as 


Here we show that the phase estimation algorithm can 
be used to implement C satisfying Eq. (HU). Let the 
mainspace be in an eigenstate |A) of the operator S with 
the corresponding eigenvalue Here |A) can be |A±) or 
|Aj_,/i) for any h. We choose the workspace Hu, to be the 
Hilbert space of /x qubits. The initial state of workspace 
\a) is chosen to be the state in which all /x qubits are in 
|0) state. We call this state all-zero state and we denote 
it by |0'). Thus our initial state is |A,0') = |A)|0'). 

The operator V corresponding to the phase estimation 
algorithm is a successive application of three operators 
in the joint space, given by 

iP = (ljv®^^)(c.5")(ljv®fT)|A,0'). (22) 

The first operator Vi — In (S' W applies the Walsh- 
Hadamard transform on the workspace and leaves the 
mainspace unchanged. Walsh-Hadamard transform is 
application of Hadamard gate on all /x qubits and it trans¬ 
forms | 0 ^) state to a uniform superposition of all basis 
states \z) (z G {0,1,..., 2^ — 1}) of the workspace. We 
have 


2^ — 1 

= ^ ^ |A)|z). (23) 

z=0 

The second operator V 2 = CzS^ is a controlled applica¬ 
tion of operator on the mainspace. It applies z itera¬ 
tions of S on the mainspace if and only if the workspace 
is in the |z) state. We get 

.. 2'*-! 

^ 2 Hi|A, 0 ') = ^ ^(5^|A))|z). (24) 

2=0 

As the mainspace is in an eigenstate |A) of 5 and the 
corresponding eigenvalue is e®'’', we get 




1 


2 '*-! 

k,z—0 



2TTzk\ 
2^ ) 


\k). 


(28) 


This is the standard output state of the phase estimation 
algorithm and has been analyzed well in literature (for 
example, see Sec. 5.2.1 of [ll|). To get an idea, note 
that we have 


.. 2'*-! 

= {exp[X - {2TTk/2f^)]y , (29) 

2=0 


which is the sum of a geometric series, and after little 
calculation we get 


|(fc|<^A)| 


1 8^(2^^ ^A —Trfc) 

2^ sin[(2^“iA — 7rA:)/2^] 


(30) 


What makes the state \(j>\) so useful is the fact that it 
provides a very good estimate of the eigenvalue of 
S provided the mainspace is in corresponding eigenstate 
I A). The accuracy of this estimation increases with the 
increasing value of /x, the number of qubits used for phase 
estimation algorithm. 

To elaborate it, we notice that each basis state \k) 
of the workspace corresponds to an integer number k G 
{ 0 , 1 ,..., 2 ^ — 1 } which is the numerical value of the bi¬ 
nary number encoded by the basis state. Let Prob(k) 
denote the probability of getting the basis state |fc) af¬ 
ter measuring the state As shown in Sec. 5.2.1 

of [HI, the probability distribution Prob(k) has a very 
sharp peak at k = kx, where 


, nearest A , 

kx ~ 2'"—(mod2^) , (31) 

Ztt 


2^-1 

^27’i|A, 0') = ^ ^ e*^^|A)|z). (25) 

2=0 

The third operator H 3 = ( 8 > -A'*' applies the inverse 

quantum fourier transform on the workspace but leaves 
the mainspace unchanged. The action of in a 2^- 
dimensional Hilbert space on each basis state |z) is given 

by 

= ^ E (26) 

k=0 

where k is the numerical value of the binary num¬ 
ber represented by the bit string encoded by the basis 
state |fc). For example, for |fc) = |0'), all qubits are 
in |0) state and hence k = 0. As suggested by Eq. 
(EH), before applying the workspace is in the state 
( 1 / 2 ^/^) Ez ea:p(xAz)|z). So, using (I26L we get 

H|A,0') =H 3 H 2 Hi|A, 0 ') = \X)\cl>x), (27) 


where the symbol ~ denotes that L.H.S. is the nearest 
integer to the number on R.H.S. The probability Prob(k) 
falls very fast as we go away from its peak at k = kx- 
To quantify this, let c be a positive integer. We define a 
subspace of the workspace Hw spanned by those basis 
states |fc) of Hu, whose numerical value k ranges from 
{kx — c)(mod2^) to {kx + c)(mod2^) in the increasing 
order. 

For clarity, we give an example to elaborate on the 
increasing order. Suppose 2>^ = 64, = 3, c = 

5. Then {kx — c)(mod2'^) = (3 — 5)mod64 = 62 
and {kx + c)(mod2^) = (3 -I- 5)mod64 = 8. By in¬ 
creasing order, we mean that the basis states \k) of 
the subspace y^ for this example encode the numbers 
(62,63, 64,1, 2,..., 8}, not the other way around like 
(62,61, 60,..., 9,8} which is the decreasing order. On 
the other hand, if kx = 62 in this example, then {kx — 
c)(mod2^) = (62 —5)mod64 = 57 and (fcA+ c)(mod2^) = 
(62 -I- 5)mod64 = 3 and the basis states |fc) of encode 
the numbers (57, 58,..., 64,1, 2, 3} not (57, 56,..., 4, 3}. 
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We now refer to Eq. (5.34) of Sec. 5.2.1 of [Ill- 
can be translated in our notation as 

Prob(k e 3^^) = > 1 - 1 ^ ■ 

(32) 

where |( 3 ^^|(/)a)| denotes the magnitude of amplitude of 
\(j)\) state in subspace and Prob (k S 3^^) denotes the 
probability of getting a basis state of 3 ^^ after measuring 
\(j)\) state. 

C. Distinguishing eigenstates 

To understand, how the phase estimation algorithm 
helps us to distinguish |A±) from |A_l), we refer to the 
discussion after Eq. dS]) at the beginning of this section. 
In our analysis of general quantum search algorithm, the 
assumption |A±| ^ 0min is crucial and we can write 

|A±| < (l-2^)0„i„ , <5<1, 

|A±| > . (33) 

Here the last inequality is due to the fact that A± are the 
only eigenvalues of S in the interval [— 6 *min; We 

now define 

^ nearest , ^ ^ , 

9=il- 6)e^in , k ^ 2'^0/27r . (34) 

With this definition and Eq. dMl), it is easy to check the 
following equations. 

0 ^ {-9,9} |A± - d| > (50min , 

9 e [-9,9} |Ai -9\> 59^,^ (35) 


Suppose that the mainspace is in |A±) eigenstate of S. 
Then Eqs. (15^ . (1571) and (1551) imply that the output state 
of workspace after phase estimation algorithm | 0 a±) is 
such that 

|(df|(^A±)|" > |(3 ^I±|<(>a±)|" > 1 - , (39) 

^ ^min 

where we have assumed 1. In above 

equations, |(<T|^a±)I is the component of |^!>a±) in df- 
subspace. Using the fact \{X\(j)x±)\^ + |(■T■*■|^A±)P = 1 
and above equation, we get 

Similar arguments show that if the mainspace is in |Aj_) 
state, we get the inequality 

In above two equations, we have not explicitly referred to 
the state of mainspace as that is implied in the definition 
of the workspace state (|(()a±), I'('ax)) given by Eq. (1271) . 
Assume that we know 0min- Then we also know 9 and k 
and thus we can easily distinguish the subspace X from 
its complementary subspace X^. 

We choose Z = X and = X^ as then the criteria of 
easy distinguishability is satisfied. Then Eq. (1T51) implies 
that the operator V corresponding to phase estimation 
algorithm can be chosen as the operator C with \a) = |0'). 
Using Eqs. ( 001 ) and m, the corresponding e is given 
by 


Let A be the subspace of Hw spanned by those 
basis states \k) whose numerical value k ranges from 
—A:(mod2'^) to fc(mod2^) in the increasing order. Let 
X^ be the complementary subspace of X. Eqs. dSU, 
(|551) and (155|) imply that for all \k) ^ X, k differs from 
^Ai by at least the numerical value of 


7 = 


27r 


(36) 


By definition, such \ k) are also the basis states of (d^Ai)"*"’ 
which is the complementary subspace of As \k) ^ 

X |fc) ^ 3^I±, we have \k) G y1^ => \k) G X and 
hence is a subspace of X. Using similar arguments, 
we can show that yjj_ is a subspace of X^. We write 
these facts as 


C df , 3^^^ C (37) 

The value of 7 depends upon <5 but the choice of S is 
arbitrary as long as the condition 25 <C 1 is satisfied. We 
choose 6 = 27r/128 to get 

6 = 27r/128 ^ 7 = . (38) 


Thus the approximation error decreases and hence the ac¬ 
curacy increases as the number of ancilla qubits used for 
the phase estimation algorithm is increased. As we need 
0{B) approximations, each approximation must have an 
error e = 0(l/i?). Then Eq. (1421) implies that 



Each application of the operator V or uses 2^ applica¬ 
tions of the operator S to implement the operator CzU^ 
as defined in Eq. (1551) . Hence the operator TZ (or ^aJ 
uses 2^^+^ applications of S and equal number of oracle 
queries as each S requires an oracle query to implement 
It- We need 0{B) applications of to transform jic) 
state to the target state \t). With the value of 2^ given 
by Eq. (l43l) . the total number of queries used in our 
post-processing is 


O 



B^ 


= 0 


B^ 



(44) 
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The number of ancilla qubits required by our algorithm 
is the logarithm of Eq. (H51) . i.e. 

= 2 log 2 B - log 2 Omin + b, (45) 

where 6 = 0 ( 1 ) is a small constant. 

D. Improving the query complexity 

We show that the query complexity of our post¬ 
processing algorithm can be improved by a factor of 
0 (i?^/InB) from 0{B^/Omin) to O{BliiB/0mm)- We 
point out that in case of i? 1 , we don’t need any 
post-processing as then w-state has a significant overlap 
0(1/5) with the target state. So for the cases of interest, 
5^1 and hence improving the query complexity by a 
factor of 0(5^/In 5) is important. To get this improve¬ 
ment, we don’t choose /i, the number of qubits for phase 
estimation algorithm, as given by Eq. (gSl). Rather, we 
choose it to be 


^l = - log 2 Omin + 16. (46) 

Putting this choice in Eq. gives us e = 0(2“®) which 
is a constant. If we stick to the phase estimation algo¬ 
rithm then the only way to further reduce this error to its 
tolerable value 0(1/5) is to add more qubits for phase es¬ 
timation. Eq. implies that adding 2 more qubits for 
phase estimation will halve the value of e. Hence we add 
2 log 2 5 qubits as in Eq. (HSl) to make e = 0(1/5). Also, 
adding one qubit doubles the query complexity of phase 
estimation and hence the query complexity increases by 
a factor of 5^ just to reduce e from 0(2“®) to 0(1/5). 
We show that once we achieve e = 0(2“®) using phase 
estimation, adding more qubits for phase estimation is 
not an efficient way to further reduce e. Rather, a bet¬ 
ter way is available which uses extra qubits in a more 
efficient way. 

We note that we are interested in only two mutually 
complementary and distinguishable subspaces X and X^ 
of the workspace. Let measurement of |(/a) state output a 
state within these subspaces X and X-^ with a probability 
of ProbA {X) and 1 — ProbA {X) respectively. With the 
choice (HSl) of fi in Eqs. (HOI) and dm), we have 

ProbAi = ProbAi (A) = = 0(2“!°) 

ProbAi (X) = ProbAj. (X^) = 1 - 0(2“!°). (47) 

Suppose we have a single qubit in a state |/3a) which 
simulates this measurement statistics which means that 
measuring |/3a) state gives the output 1 with the proba¬ 
bility ProbA (X) and 0 with the probability 1—ProbA {X). 
We write it as 

\/3f) = cos/3a| 0) ± *sin/3A|l), sin/3A = | VProbA iX)\. 

(48) 

We could have chosen any phase factor for the |1) state 
in above equation but the choice ± is sufficient for our 


purpose. The above single qubit state is a very important 
resource to improve our approximation. To understand 
this, suppose we have v identical copies of such a qubit. 
Let 52 *' denote the joint Hilbert space of all i' qubits. 
The state of v identical copies of |/3^) is given by 

I/3a ,^) = (cos/3a| 0) ±zsin/3A|l))®'' • (49) 

We now partition the Hilbert space 'H 21 ' in two mutually 
complementary subspaces: W and defined as 

(v/l)-! 

w = U |0)"'|1)"“"' , 

V 

W-^ = (J | 0 )‘''| 1 )'^“‘'' . (50) 

v'=v/2 


Here 10)“" denotes that subspace of 'H 2 '' in which 

exactly v' qubits are in | 0 ) state while remaining v — v' 
qubits are in 11 ) state. 

We now find |(yV|/3^, z^)|, the magnitude of component 
of the joint qubit state \f5'^,v) in the subspace W. Its 
square is the probability of getting a state within W sub¬ 
space, i.e. of getting at least half of v qubits in | 1 ) state 
after measuring |/3^,i/). Let Prob(i/') denote the prob¬ 
ability of getting v' qubits in | 1 ) state. As the state of 
V qubits is a tensor product of single qubit states, the 
probability distributions of measurement outcomes of all 
qubits are independent of each other. Hence Prob(i^') 
follows the binomial distribution. Eq. (1481) implies that 
the probability of getting a single qubit in | 1 ) state is 
sin^ j3\ = ProbA {X). Then the binomial distribution im¬ 
plies that Prob(j^') has a sharp peak at v' = i/x ProbA (X) 
and decays exponentially with increasing v' as we go 
away from this sharp peak. Eq. (1471) implies that this 
sharp peak is at i^ProbA (X) = 1 / [l — 0(2“®°)] fn v for 
A = A± and i^ProbA (X) = z/0(2“®°) for A = A_l. Hence 

V 

A = A± => ^ Prob(j^') = 1 — 0 ( 6 ““^) (51) 

v'=v/2 


and similarly A = Aj_ implies that the quantity on R.H.S. 
is o(e“^). Here o{x) <C a: is the small-o notation. As 
mentioned earlier, this quantity is square of |(W|,5^, i')\. 
Hence 


A = A± |(W|/3±±,;.)| = l-o(e“'^), 

A = A^ = (52) 

By choosing 

v = b\nB, (53) 


we get 

A = A± ^ |(W|/3±±,z^)| = l-0(1/5), 
A = A^ |(W|/3±^,z^)| =0(1/5), 


( 54 ) 
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Thus there exists a unitary operator C satisfying the 
properties of Eq. (HU) with e = o{\/B) if we choose the 
subspaces to be Z = W and If we use this 

operator to approximate , we get a negligible error in 
our algorithm. It requires preparation of v qubits in the 
state 1/3^) given by (l4^ . The question arises: how do we 
prepare such a qubit? As it simulates the measurement 
statistics of the workspace state |^a), naturally the pro¬ 
cess of getting 1/3^) state has to depend upon the process 
of getting \cj)x) state which is the phase estimation algo¬ 
rithm. Consider the following amplitude amplification 
operator acting on the workspace 

A{,cj^x,X) = -I^Jx (55) 

where is the selective phase inversion of |((ia) state 
and /a: is the selective phase inversion of X subspace. 
Ix is implemented using our knowledge of the subspace 
X. To implement we note that this is equivalent to 
the selective phase inversion of |A)| 0 a) as the mainspace 
is in I A) state. Now, Eq. (l27l) gives us 

iP|A,0') = |A)|</.A), (56) 

which implies that the selective phase inversion of |A)|(^a) 
is equal to where Ia.o' is the selective phase 

inversion of |A)|0') state and it is equivalent to imple- 
mentable /q', the selective phase inversion of | 0 ') as the 
mainspace is in |A) state. As discussed earlier, both V 
and pi requires 2^ queries for implementation. Thus the 
operator A{(j)\,X) requires 2^^+^ queries. 

The eigenspectrum of amplitude amplification opera¬ 
tor has been analysed in detail in Section 2 of Q . It has 
two non-trivial eigenstates given by 

|A±) = (1/^2) (|Aa) ± \X^)) , (57) 

where \X\) (iTj;^)) is the projection of \cj)\) state on X 
{X^) subspace. The corresponding eigenvalues are given 

by 

Wa = Sin“^ | (((lA |-Ta) | . (58) 

Comparing above equation with Eq. (1481) implies that 

|wa| = I^aI (59) 

Also, it was shown in Q that the state \cj)\) is a linear 
combination of eigenstates \X'^). Precisely 

|,^a) = (I/V 2 ) \X+) - e-*-- \X-)) (60) 

Suppose our workspace is in an eigenstate \X^) of 
A{4>x,^) which we now denote simply by A for con¬ 
venience. To get a qubit in |/3^) state, we attach a qubit 
initially in the state |0) to our workspace. We apply a 
Hadamard gate on it to get the joint state of workspace 
and qubit as 

|A+)(^(| 0 ) + | 1 ))) . (61) 


We now apply a controlled-operation which applies the 
operator A on the workspace if and only if the qubit is 
in |1) state. Then we get the following joint state 

^(|A+)|0)-byl|A+)|I)) . (62) 

As \X^) is an eigenstate of A with the eigenvalue 
the above state is 

^{\X+m + A^-^\X+)\l)) . (63) 

which is equivalent to 

|A+)(^^(|0) + e*2-Ml))) . (64) 

We now apply a Hadamard gate on the qubit and the 
above state becomes, upto a global phase, 

lA’;^^) (coswa| 0) — zsinwAll)) • (65) 

The qubit state is same as |/3^) as |wa| = I/SaI due to Eq. 

(EH). 

As the workspace state remains unchanged in the 
above operation (it remains in the eigenstate jd?^) of 
A), we can attach arbitrary, say l', number of qubits and 
repeat the above operation i' times using 1 / applications 
of A to get the state 

l-^A^) (I/3 a“))®'^ ■ (66 ) 

Similar arguments can show that if the workspace is in 
other eigenstate \X^) of A^ then we get the following 
state 

I^a") (1/31))®'^ • (67) 

But in our case, the workspace is in \(p\) state which is 
a linear superposition given by Eq. dsni)- Hence doing 
above operations will give us the state 

^ (l^A^) (I/3 a-))®^ + (l/3|))®0 (68) 

where r is a phase factor irrelevant for our purpose. This 
linear combination of |/3^) states is as good as or 
1/3^) for our purpose as it doesn’t change the measure¬ 
ment statistics of getting O’s or I’s. 

Now we are in a position to determine the query and 
space complexity of the operator TZ. To implement 77., 
we need 

(i) One copy of |^a) state in the workspace. It takes 
2 ^ = O(l/0min) queries and /r = — log 2 6*niin + 16 qubits. 

(ii) We then need = 5 In 73 applications of the operator 
A{(j)\,X). As each A takes 2^^^ queries, this requires 
0(ln732'^) = O(lni3/0min) queries and 5 In 73 qubits. 

We need 0(73) applications of 77 for our algorithm to 
evolve \w) state to the target state \t). We add 0{B/a) 


queries, required to get the [w) state, to get the total 
query complexity of our algorithm as 


O 



B\^B 

^min 


(69) 


The number of ancilla qubits required by our algorithm 
is 


^ = - log 2 6»min + 51ni? + 16. (70) 

This is almost same as given by Eq. (1451) which uses only 
phase estimation algorithm to achieve the tolerable error 
in approximation. This concludes our discussion on the 
approximation algorithm for selective phase inversion. 


IV. DISCUSSION AND CONCLUSION 

We have shown in this paper that there exists a better 
way to transform the output state |w) of general quan¬ 
tum search algorithms to the target state \t) which en¬ 
codes the solution of search problems. Many important 
variants of the Grover’s search algorithm can be shown to 
be special cases of the general quantum search algorithm 
as discussed in Section III of Not just that, the 

general algorithm can also be used to design quantum 
search algorithms specifically suited for certain compu- 
tational problems as illustrated in a recent paper [13]. 
Hence we believe that the post-processing algorithm pre¬ 
sented in this paper can find important applications in 
search problems. 

In Section III of this paper, we have also presented an 
algorithm to efficiently approximate the selective phase 
inversion of unknown eigenstates of any operator. We 
believe that our approximation algorithm is important 
in its own context and can find useful applications. For 
example, in a recent paper 0, same algorithm is used 
in the context of simulating the quantum adiabatic evo¬ 
lution [Ts] on a digital quantum computer. This simula¬ 
tion is faster than earlier work [l^. We can also use our 
approximation algorithm to improve the time complex¬ 
ity of adiabatic quantum computation with one dimen¬ 
sional projector Hamiltonian [^. There, the final state 
of quantum adiabatic evolution has an overlap of 1/C 
with the desired target state, where C depends upon the 
eigenspectrum of the initial Hamiltonian in a similar way 
as B depends upon the eigenspectrum of Dg in our al¬ 
gorithm. This final state, denoted as \E_^^+) in [a, 
can be transformed to the target state using quantum 
amplitude amplification provided we can implement the 
selective phase inversion of |i?_ ^+). Our approximation 


algorithm can be used for this implementation as |i?_ ^+) 
is the ground state of a known Hamiltonian and it can be 
reliably distinguished from excited states using the phase 
estimation algorithm. 

Referring to the discussion after Eq. m, it appears 
that we need to know 0niin in advance to use our approx¬ 
imation algorithm. This knowlegde is required to distin¬ 
guish I A±) eigenstates of S from other eigenstates by esti¬ 
mating their eigenvalues as |A±| <C 0min and |Aj_) > 0min- 
Hence 0niin needs to be known to do any selective trans¬ 
formation on |A±). Suppose we don’t know 6*niin and by 
mistake, we choose it to be such that apart from A±, 
there also exists non-A± eigenphases of S in the interval 
{—fi'min, ^*min}- In this case, our algorithm applies a se¬ 
lective phase inversion on not just the |A±) eigenstates 
but also some of the non-|A±) eigenstates. Similarly, if 
we choose ^min to be very small so that no eigenphase 
of S lies in the interval {—6)min, 6*niin} then our algorithm 
fails to do any selective transformation and implements 
a close to identity operation. In both these cases, our 
algorithm doesnot have any success gaurantee. 

In case of no prior knowlegde of Omin, we need to run 
multiple, say r, rounds of our algorithm. First, we begin 
by choosing 0min to be say If it is larger than 

actual 0min, then the algorithm may fail and in that case, 
we apply the algorithm again by choosing a decreased 
value (1 — We keep on doing so and at (r)**' 

application of the algorithm, the chosen value of 0niin is 
(1 — :^)’'2“^°. For r — 0(— In^min), the chosen value will 
be sufficiently close to the actual value ^min to provide a 
successful quantum search algorithm. 


For the example of two-dimensional spatial search 
o = ^min = ^/VN and B = 0{\ogN) and hence 
our algorithm is just another alternative algorithm with 
a slightly increased time complexity by a factor of 
O(loglogA^). We point out that Ambainis et.al. pre¬ 
sented an 0{^/N log N) algorithm for two-dimensional 
spatial search based on classical post-processing [l8| . Ba¬ 
sically they had shown that in case of 2d spatial search, 
w-state jw) is such that with high probability, mea¬ 
suring it gives a basis state lie') within the distance 
of X local steps from the target state 

\t). Hence, |t) can be found by classically searching a 
X 7V^/^) lattice centred at jw') which takes 0{N) 
time steps. Our algorithm offers an alternative based 
on quantum post-processing and may be more suitable if 
two-dimensional spatial search is used as a quantum sub¬ 
routine for other algorithms. However, we point out that 
our algorithm mainly offers advantages only for general 
quantum search algorithms where 9min/ log B ^ a which 
is a relaxed criteria. 
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